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Electric and magnetic fields of fractal distribution of charged particles are consid- 
ered. The fractional integrals are used to describe fractal distribution. The fractional 
integrals are considered as approximations of integrals on fractals. Using the frac- 
tional generalization of integral Maxwell equation, the simple examples of the fields of 
homogeneous fractal distribution are considered. The electric dipole and quadrupole 
moments for fractal distribution are derived. 



I. INTRODUCTION 

Derivatives and integrals of fractional order have found many applications in recent 
studies in physics. The interest in fractional analysis has been growing continually in recent 



years. Fractiona, 
mechanics 



analysis has numerous applications: kinetic 



,2. 



theories 12, 3, 41; statistical 



Bfl Q; dynamics in complex media fl U 

[ill \\\ ; electromagnetic theory 
I3I I4I 0,^6 1 and many others. The new type of problem has increased rapidly in areas 
in which the fractal features of a process or the medium impose the necessity of using 
nontraditional tools in "regular" smooth physical equations. In order to use fractional 
derivatives and fractional integrals for fractal distribution, we must use some continuous 
.ned.un, .node, fl . We propose to describe the fracta, d.stributiou b y a fraction*, continuous 
medium [jj, where all characteristics and fields are defined everywhere in the volume but 
they follow some generalized equations which are derived by using fractional integrals. In 
many problems the real fractal structure of medium can be disregarded and the fractal 
distribution can be replaced by some fractional continuous mathematical model. Smoothing 

of microscopic characteristics over the physically infinitesimal volume transforms the initial 

n 

fractal distribution into fractional continuous model |9( that uses the fractional integrals. The 
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order of fractional integral is equal to the fractal dimension of distribution. The fractional 
integrals allow us to take into account the fractality of the distribution. Fractional integrals 
are considered as approximations of integrals on fractals In Ref. Q], the authors 

proved that integrals on net of fractals can be approximated by fractional integrals. In 
Ref. |§], we proved that fractional integrals can be considered as integrals over the space 
with fractional dimension up to numerical factor. In order to prove, we use the formulas of 
dimensional regularizations [3]. 



;> articles. 
9. 
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We can consider electric and magnetic fields of fractal distribution of chargec 
Fractal distribution can be described by fractional continuous medium model 
In the general case, the fractal distribution of particles cannot be considered as continuous 
medium. There are points and domains that have no particles. In Ref. |9j , we suggest to con- 
sider the fractal distributions as special (fractional) continuous media. We use the procedure 
of replacement of the distribution with fractal mass dimension by some continuous model 
that uses fractional integrals. This procedure is a fractional generalization of Christensen 
approach 2l3|. Suggested procedure leads to the fractional integration and differentiation 
to describe fractal distribution. In this paper, we consider the electric and magnetic fields 
of fractal distribution of charged particles. In Sec. II, the densities of electric charge and 
current for fractal distribution are considered. In Sec. Ill and Sec. IV, we consider the 
simple examples of the fields of homogeneous fractal distribution. In Sec. V, we consider 
the fractional generalization of integral Maxwell equation. In Sec. VI, the examples of elec- 
tric dipole and quadrupole moments for fractal distribution are considered. Finally, a short 
conclusion is given in Sec. VII. 



II. ELECTRIC CHARGE AND CURRENT DENSITIES 
A. Electric charge for fractal distribution 

Let us consider a fractal distribution of charged particles. For example, we can assume 
that charged particles with a constant density are distributed over the fractal. In this case, 
the number of particles N enclosed in a volume with characteristic size R satisfies the scaling 
law N(R) ~ R D , whereas for a regular n-dimensional Euclidean object we have N(R) ~ R n . 

For charged particles with number density n(r,t), we have that the charge density can 
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be defined by 

p(r,t) = qn(r,t), 

where q is the charge of a particle (for electron, q = — e). The total charge of region W is 
then given by the integral 

Q(W)= [ P (r,t)dV 3 , 
Jw 

or Q(W) = qN(W), where N(W) is a number of particles in the region W. The fractional 
generalization of this equation can be written in the following form 



Q(W)= / p(r,t)dV D , 
Jw 

where D is a fractal dimension of the distribution, and dVr, is an element of D-dimensional 
volume such that 

dV D = c 3 (D,r)dV 3 . (1) 

For the Riesz definition of the fractional integral jjj, the function c 3 (D, r) is defined by 
the relation 

23-^X3/2) 

C3(Ar) " T(D/2) M • (2) 
The initial points of the fractional integral are set to zero. The numerical factor in Eq. (J2J) 

has this form in order to derive usual integral in the limit D — > (3 — 0). Note that the 

usual numerical factor %\D) = T(l/2)/[2 D 7c 3 / 2 T(D/2)], which is used in Ref. Q leads to 

7 - 1 (3 - 0) = r(l/2)/[2 3 vr 3 / 2 r(3/2)] in the limit D -»• (3 - 0). 

For the Riemann-Liouville fractional integral [1], the function c 3 (D, r) is defined by 

Here we use Cartesian's coordinates x, y, and z. In order to have the usual dimensions of 
the physical values, we can use vector r, and coordinates x, y '. ^ cis dimensionless values. 

Note that the interpretation of fractional integration is connected with fractional di- 
mension [^, 0]. This interpretation follows from the well-known formulas for dimensional 
regularizations . The fractional integral can be considered as an integral in the fractional 
dimension space up to the numerical factor T(D / 2) /[27r D / 2 T(D)}. 

If we consider the ball region W = {r : |r| < R}, and the spherically symmetric 
distribution of charged particles (n(r,t) = n(r)), then we have 

N(R) = 4tt W > ^ n(r)r D ^dr, Q(R) = qN(R). 



4 



For the homogeneous (n(r) = n ) fractal distribution, we get 

Fractal distribution of particles is called a homogeneous fractal distribution if the power law 
N(R) ~ R D does not depend on the translation of the region. The homogeneity property 
of the distribution can be formulated in the following form: For all regions, W and W 
such that the volumes are equal, V^W 7 ) = V(W), we have that the numbers of particles 
in these regions are equal too, N(W) = N(W). Note that the wide class of fractal media 
satisfies the homogeneous property. In Ref. p|, the continuous medium model for the fractal 
distribution of particles was suggested. Note that the fractality and homogeneity properties 
can be realized in the following forms: 

(1) Homogeneity: The local number density of homogeneous fractal distribution is a trans- 
lation invariant value that has the form n(r) = n = const. 

(2) Fractality: The number of particles in the ball region W obeys a power law relation 
N D {W) ~ R D , where D < 3, R is the radius of the ball. 



B. Electric current of fractal distribution 



For charged particles with number density n(r,t) flowing with velocity u = u(r, t), the 
resulting density current J(r, t) is given by 

J(r,t) = gn(r,t)u, 

where q is the charge of a particle (for electron, q = —e). 

The electric current is defined as the flux of electric charge. Measuring the field J(r, t) 
passing through a surface S = dW gives the current (flux of charge) 

I(S) = $j(S)= [ (J,dS 2 ), 

where J = J(r, t) is the current field vector, (IS2 = dS^n is a differential unit of area pointing 
perpendicular to the surface S, and the vector n = n^ek is a vector of normal. The fractional 
generalization of this equation for the fractal distribution can be written in the following 
form: 

I{S)= [ (J(r,t),rfS d ), 
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where we use 

22-d 

dS d = c 2 (d, r)dS 2 , c 2 (d, r) = |if~ 2 . (4) 

Note that c 2 (2, r) = 1 for d — 2. In general, the boundary dW has the dimension d. In the 
general case, the dimension d is not equal to 2 and is not equal to (D — 1). 



C. Charge conservation for fractal distribution 

The electric charge has a fundamental property established by numerous experiments: 
the change of the quantity of charge inside a region W bounded by the surface S = dW is 
always equal to the flux of charge through this surface. This is known as the law of charge 
conservation. If we denote by J(r, t) the electric current density, then charge conservation 
is written 

dQ(W) 

or, in the form 



d 
~dt 

In particular, when the surface S = dW is fixed, we can write 



/ p{T,t)dV D = - I (J(r,t),dS d ). (5) 
Jw Jaw 



Using the fractional generalization of the mathematical Gauss's theorem (see Appendix), 
we have 

/ (3(r,t),dS d )= [ c^(D,r)^-(c 2 (d,r)J k (r,t))dV D . (7) 
Jaw Jw dx k V / 

Substituting the right hand sides of Eqs. (JE1) and ((Z|) in Eq. 0, we find the law of charge 
conservation in differential form 

c 3 (Ar)^|^ + ^-(c 2 (rf,r)J fc (r,t)) =0. 

This equation can be considered as a continuity equation for fractal distribution of particles 

0. 
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III. ELECTRIC FIELD OF FRACTAL DISTRIBUTION 

A. Electric field and Coulomb's law 

For a point charge Q at position r' (i.e., an electric monopole), the electric field at a point 
r is defined in MKS by 

E Q r " r ' 



Attsq |r — r'| 3 ' 

where Eq is a fundamental constant called the permittivity of free space. 

For a continuous stationary distribution p(r') of charge, the electric field E at a point r 
is given in MKS by 

where Eq is the permittivity of free space. For Cartesian's coordinates dV% = dx'dy'dz' . 

The fractional generalization of Eq. (jHJ) for a fractal distribution of charge is given by 
the equation 

where dV' D = c${D, r')dV^. Equation (jHJ) can be considered as Coulomb's law written for a 
fractal stationary distribution of electric charges. 

Measuring the electric field passing through a surface S = dW gives the electric flux 

$ E (S)= [ (E,dS 2 ), 
Js 

where E is the electric field vector, and dS2 is a differential unit of area pointing perpendic- 
ular to the surface S. 

B. Gauss's law for fractal distribution 

Gauss's law tells us that the total flux <&e{S) of the electric field E through a closed 
surface S = dW is proportional to the total electric charge QiW) inside the surface: 

§ E (dW) = -Q{W). (10) 

The electric flux out of any closed surface is proportional to the total charge enclosed within 
the surface. 
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For the fractal distribution, Gauss's law states 

/ (E,dS 2 ) = - / p{r,t)dV D (11) 
Js £ o Jw 

in MKS, where E = E(r, t) is the electric field, and p(r,t) is the charge density, dVo = 
c 3 (D, r)dV 3 , and Eq is the permittivity of free space. 

Gauss's law by itself can be used to find the electric field of a point charge at rest, and the 
principle of superposition can then be used to find the electric field of an arbitrary fractal 
charge distribution. 

If we consider the stationary spherically symmetric fractal distribution p(r, t) = p(r), and 
the ball region W = {r : |r| < R}, then we have 

r R 

Q(W)=4n p(r)c 3 (D,r)r 2 dr, 
Jo 

where c 3 (D,r) is defined in Eq. (J2J), i.e., 

Q{w) = a * ^t°d% 2) [p^y D ' ldr - ( 12 ) 

Using the sphere S' = {r: |r| = _R}asa surface S = dW, we get 

^ E (dW) =AttR 2 E(R). (13) 

Substituting Eqs. (fT^j) and (|T3*j) in Gauss's law (|IDj). we get the equation for electric field. 
As a result, Gauss's law for fractal distribution with spherical symmetry leads us to the 
equation for electric field 

E(R) = 23 ~!£ (3/2 l I* p{r)r D - l dr. 

For example, the electric field of homogeneous (p(r) = p) spherically symmetric fractal 
distribution is defined by 

771^ _ ^ ^(3/2) D _2 tdD-2 

E(R) ~ P e DT(D/2) R ~ R ■ 

IV. MAGNETIC FIELD OF FRACTAL DISTRIBUTION 

A. Magnetic field and Biot-Savart law 

The Biot-Savart law relates magnetic fields to the currents which are their sources. In 
a similar manner, Coulomb's law relates electric fields to the point charges which are their 
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sources. Finding the magnetic field resulting from a fractal current distribution involves the 
vector product and is inherently a fractional calculus problem when the distance from the 
current to the field point is continuously changing. 

For a continuous distribution the Biot-Savart law in MKS has the form 

where [ , ] is a vector product, J is the current density, /x is the permeability of free space. 
The fractional generalization of Eq. (|14|) for a fractal distribution in MKS has the form 

W-ZjJ^Sfa ( 15 ) 

This equation can be considered as Biot-Savart law written for a steady current with fractal 
distribution of electric charges. The Biot-Savart law (|THj) can be used to find the magnetic 
field produced by any fractal distribution of steady currents. 



B. Ampere's law for fractal distribution 

The magnetic field in space around an electric current is proportional to the electric 
current which serves as its source, just as the electric field in space is proportional to the 
charge which serves as its source. In the case of static electric field, the line integral of the 
magnetic field around a closed loop is proportional to the electric current flowing through the 
loop. The Ampere's law is equivalent to the steady state of the integral Maxwell equation 
in free space, and relates the spatially varying magnetic field B(r) to the current density 

J(r) ' n 

Note that, as mentioned by Lutzen in his article |21lj . Liouville, who was one of pioneers 
in development of fractional calculus, was inspired by the problem of fundamental force law 
in Ampbre's electrodynamics and used fractional differential equation in that problem. 

Let be a closed path around a current. Ampere's law states that the line integral of the 
magnetic field B along the closed path L is given in MKS by 

£{B,dl) = f i I{S) i 

where dl is the differential length element, and /xq is the permeability of free space. For the 
fractal distribution of charged particles, we use 

I(S)= [ (J,dS d ), 
Js 
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where dSj = c 2 (d, r)dS 2 . 

If we consider the cylindrically symmetric fractal distribution, we have 

I(S) = 2tt / J(r)c 2 (d,r)rdr, 
Jo 

where c 2 (d, r) is defined in Eq. (jlj, i.e., 

/ < s) = fc r|J)/' ij(r)r " 1 *' 

Using the circle L = dW = {r : |r| = R}, we get 

(B,dl) = 27T.R 5(7?). 

L 

As a result, Ampere's law for fractal distribution with cylindrical symmetry leads us to the 
equation for magnetic field 



For example, the magnetic field B(r) of homogeneous (J(r) = Jo) fractal distribution is 
defined by 

V. FRACTIONAL INTEGRAL MAXWELL EQUATIONS 

The Maxwell equations are the set of fundamental equations governing electromagnetism 
(i.e., the behavior of electric and magnetic fields). The equations that can be expressed in 
integral form are known as Gauss's law, Faraday's law, the absence of magnetic monopoles, 
and Ampere's law with displacement current. In MKS, these become 

/(E,dS 2 ) = - / pdV D , 
Js £ o Jw 

jf(E,dl 1 ) = ~y , (B,dS a ) > 
^(B,dS 2 )=0, 
jf (B, dh) = fi jT(J, dS d ) + 6 0f i -^- t J(E, dS 2 ). 
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Let us consider the special case such that the fields are defined on fractal 221 only. 
The hydrodynamic and thermodynamics fields can be defined in the fractal media % [ll| . 
Suppose that the electromagnetic field can be defined on fractal as an approximation of 
some real case with fractal medium. If the electric field E(r) and magnetic field B(r) can 
be defined on fractal and does not exist outside of fractal in Eucledian space E 3 , then we 
must use the fractional generalization of the integral Maxwell equations in the form 

[E,dS d ) = - [ P dV D , 
e o Jw 

'^,dL r ) = ~ I (B,dS d ), 



[B,dS d ) = 0, 

j{B,d\ 1 ) = p J {J,dS d )+e p -^ j (E,dS d ). 

These fractional integral equations have unusual properties. Note that fractional integrals 
are considered as an approximation of integrals on fractals jl?! l"!^. 

Using the fractional generalization of Stokes's and Gauss's theorems (see Appendix), we 
can rewrite the fractional integral Maxwell equations in the form 

/ c^(D,r)div(c 2 (d,T)B)dV D = — [ pdV D , 
Jw e o Jw 

%\D, r)div{c 2 {d, r)B)dV d = 0, 

J cj x (ti, r)(curl(c 1 (7, r)B), dS d ) = fi J (J, dS d ) + £ ol^o^ J (E,dS rf ), 
As a result, we have the following differential Maxwell equations: 

div(c 2 (d, t)E) = —c 3 (D,r)p, 

cwr/(ci(7,r)E) = -c 2 (d,r)^B, 

div(c2{d, r)Bj = 0, 

/ \ 9E 

curl [ci (7, r)B J = p c 2 (d,r)J + e p c 2 (d } r)—. 
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Note that the law of absence of magnetic monopoles for the fractal leads us to the equation 
div(c 2 (d, r)B) = 0. This equation can be rewritten in the form 



In the general case (d ^ 2), the vector grad (c 2 (d, r)) is not equal to zero and the magnetic 
field satisfies div~Q ^ 0. If d — 2, we have div (F) ^ only for nonsolenoidal field F. 
Therefore the magnetic field on the fractal is similar to the nonsolenoidal field. As a result, 
the magnetic field on fractal can be considered as a field with some "fractional magnetic 
monopole" q m ~ (B, Vc 2 ). 

VI. MULTIPOLE MOMENTS FOR FRACTAL DISTRIBUTION 
A. Electric multipole expansion 

A multipole expansion is a series expansion of the effect produced by a given system in 
terms of an expansion parameter which becomes small as the distance away from the system 
increases. Therefore, the leading one of the terms in a multipole expansion are generally 
the strongest. The first-order behavior of the system at large distances can therefore be 
obtained from the first terms of this series, which is generally much easier to compute than 
the general solution. Multipole expansions are most commonly used in problems involving 
the gravitational field of mass aggregations, the electric and magnetic fields of charge and 
current distributions, and the propagation of electromagnetic waves. 

To compute one particular case of a multipole expansion, let R = X^k be the vector from 
a fixed reference point to the observation point, r = Xk^k be the vector from the reference 
point to a point in the body, and d = R — r be the vector from a point in the body to the 
observation point. The law of cosines then yields 



div B 



(B, gradc 2 (d, r)). 




where d = |d|, and cos 9 = (r, R)/(r\R), so 




Now define e = r/R, and x = cos 8, then 
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/ \ — 

But (1 — 2ex + e 2 ) is the generating function for Legendre polynomials P n {x) as 
ows: 

— 1/2 °° 

(l-2ex + e 2 ) =^ e n P n (x) 

n=0 



so, we have the equation 



1 1 v^/r 



oo 



n=0 

Any physical potential that obeys a (1/cf) law can therefore be expressed as a multipole 
expansion 

oo If 

The n = term of this expansion, called the monopole term, can be pulled out by noting 
that Pq(x) = 1, so 

U = i5 1 " (r) * + i £ fl^ £ ^.(cc^r)^. (17) 



The nth term 

«„ 1 1 



^ / r"i>„(cos9)p(r)<{Vi (18) 



4tt£:o 

is commonly named according to the following: n - multipole, - monopole, 1 - dipole, 2 - 
quadrupole. 

B. Electric dipole moment of fractal distribution 

An electric multipole expansion is a determination of the voltage U due to a collection 
of charges obtained by performing a multipole expansion. This corresponds to a series 
expansion of the charge density p(r) in terms of its moments, normalized by the distance 
to a point R far from the charge distribution. In MKS, the electric multipole expansion is 
given by Eq. (JTE|) : 

1 °° 1 f 

U = ^ a W»j w r^ff W)d V D , (19) 

where P n (cos#) is a Legendre polynomial and 9 is the polar angle, defined such that 

cosfl = (r,R)/(|r||R|). 

The first term arises from Po^) — 1> while all further terms vanish as a result of P n {x) 
being a polynomial in x for n > 1, giving P n (0) = for all n > 1. 
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If we have 



Q(W) = [ p(r)dV D = 0, 
Jw 



then the n = term vanishes. Set up the coordinate system so that 9 measures the angle 
from the charge-charge line with the midpoint of this line being the origin. Then the n — 1 
term is given by 

U * = 7^4 / rP 1 (cosd)p(r)dV D = 
4ire Hf J w 

For a continuous charge distribution, the electric dipole moment is given by 

vp(r)dV 3 , (20) 



4ne R 2 J w 
1 1 



' w 

where r points from positive to negative. Defining the dipole moment for the fractal distri- 
bution by the equation 

pP) = f rp(r)dV D , (21) 
Jw 

then gives 

1 (R,p (D) ) 1 p^cosa 
1 ~ A^T R~ 3 ~ A^T R 2 

where cosa = (R,p^)/(p^R), and p^ = \J {pf ] ) 2 + (p ( y D) ) 2 + {pf ] ) 2 . 

Let us consider the dipole moment for the fractal distribution by Eq. ((25), where we use 
the Riemann-Liouville fractional integral, and the function c 3 (D,r) in the form 

c a {D,r) = UW, a = D/3. (22) 

Let us consider the example of electric dipole moment for the homogeneous (p(r) = p) fractal 
distribution of electric charges in the parallelepiped region 

W = {(x; y; z) : < x < A, < y < B, < z < C}. (23) 

In this case, we have Eq. (|2ip in the form 

= ' r dx r dy r dz , v -v- - A 
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The electric charge of parallelepiped region (|2*3|) is defined by 

J w a d r d (a) 

Therefore, we have the dipole moments for fractal distribution in parallelepiped in the form 

pi D) = ^tQ(W)A, p^ = ^-Q(W)B, p ^ = ^-Q(W)C, 
a + 1 y a + 1 a + 1 

where we can use a/ (a + 1) = D/(D + 3). As a result, we get 

V?> = , (24) 

where are the dipole moments for three-dimensional homogeneous distribution. If we 
use the following limits 2 < D < 3, then we have 

2D 

0.8 < < 1. 

D + 3 ~ 



C. Electric quadrupole moment of fractal distribution 

While this is the dominant term for a dipole, there are also higher-order terms in the 
multipole expansion that become smaller as R becomes large. The electric quadrupole term 
in MKS is given by 

^ = / r 2 P 2 (cos6)p(r)dV D = 

47T£ R Jw 

J w r2 (h os26 ~l) p{r)dVD 



1 1 



Atteq R 3 

1 1 



4tT£ 2R 3 

11/'/ ■-'> 



r 2 (3cos 2 6 - l)p(r)dV D 

w 



(R, vY-r l ) P {r)dV D 



4ne 2R 3 J w \R 2 

The electric quadrupole is the third term in an electric multipole expansion, and can be 
defined in MKS by 

k,l=l 

where Eq is the permittivity of free space, R is the distance from the fractal distribution of 
charges, and Qui is the electric quadrupole moment, which is a tensor. 
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The electric quadrupole moment is defined by the equation 

Qki = / {3x k xi - r 2 5 M )p{r)dV D , 
Jw 

where x^ = x,y, or z. From this definition, it follows that 

3 

Qkl = Qlk, and y^Qfcfc = 0. 

k=l 

Therefore, we have Q zz = —Q xx — Q yy . In order to compute the values 

Qg>= / [3x 2 - (x 2 + y 2 + z 2 )]p{v)dV D = [ [2x 2 -y 2 -z 2 }p(r)dV D , 
Jw Jw 

= [ [% 2 - ^ + V 2 + z 2 )WWd = I l-x 2 + 2y 2 - z 2 )]p(r)dV D , 
Jw Jw 

Q? z ] = I [3z 2 - (x 2 + y 2 + z 2 )]p(v)dV D = [ [-x 2 -y 2 + 2z 2 )]p(r)dV D , 
Jw Jw 

we consider the following expression 

Q(a,/3, 7 )= / [ax 2 + (3y 2 + 1Z 2 )]p(r)dV D} (25) 
Jw 

where we use the Riemann-Liouville fractional integral and the function Cs(D, r) in the 
form 

C3(Ar) = ^, a = D/3. (26) 

Using Eq. (J23j) . we have 

g2 ) = 0(2, gi^ = g(-i,2,-i), gg> = g(-i,-i,2). (27) 

D. Quadrupole moment of fractal parallelepiped 

Let us consider the example of electric quadrupole moment for the homogeneous (p(r) = 
p) fractal distribution of electric charges in the parallelepiped region 

W = {{x; y; z) : < x < A, < y < B, < z < C}. (28) 

If we consider the region W in form (|28|). then we get 

(a + 2)a z L 6 (a) 
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The electric charge of this region W is 

p(ABC) a 
a 3 r 3 (a) ' 



Q(W) = p f dV D 
Jw 



Therefore, we have the following equation 

a 



Q(a, (3, 7) = — —Q{W){aA 2 + (3B 2 + 7 C 2 ), 
a + 2 

where a = D/3. If D = 3, then we have 

Q(a,(3 n ) = ^Q(W)(aA 2 + $B 2 + 1 C 2 ). 
As a result, we get electric quadrupole moments of fractal distribution in the region 



W: 

{D) _ 3D (3) 

where Q kl ! are moments for the usual homogeneous distribution [D — 3). By analogy with 



these equations, we can derive Qu?' for the case k ^ I. These electric quadrupole moments 



are 



Qi? = 77^5 off. (MO- 



(£> + 3) 2 

If we use the following limits 2 < D < 3, then we get the relations 

3D AD 2 

0.75 < < 1, 0.64 < < 1. 

D + 6 ~ ' (D + 3) 2 ~ 

E. Quadrupole moment of fractal ellipsoid 

Let us consider the example of electric quadrupole moment for the homogeneous (p(r) = 
p) fractal distribution in the ellipsoid region W: 

- + ^ + -<l. (29) 
A 2 B 2 C 2 ~ y ' 

If we consider the region W in the form (|29|) . then we get (|25|) in the form 

(3a + 2)r ,:i (a) 
where a = -D/3, and i£j(a) (i=l,2,3) are defined by 

Kx(a) = L(a + 1, a - 1, 2ix)L(a - 1, 2a + 1, tt), 
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K 2 (a) = L[a — 1, a + 1, 2tt)L(ci — 1,2a + 1, 7r), 
if 3 (a) = L(a - 1, a - 1, 2n)L(a + 1, 2a - 1, 7r). 
Here we use the following function 

21 ^ , _,_m*i ^_,_x,m « r(n/2 + l/2)r(m/2 + 1/2) 



21 r 1 

L(n,m,l) = — / dx | cos(x)| n | sin(x) 
^ Jo 



7T r(n/2 + m/2 + l) 



If D = 3, we obtain 



Q(a, /?, 7) = y + /?5 2 + 7 C 2 ), (30) 

where we use K\ — K 2 — K 3 = 47r/3. The total charge of this region W is 

Q(VP) = p / dVb = - ■ 31 

J w 3ar d (a) r(3a/2) 

If D = 3, we have the total charge 

Q(TU) = p / dU 3 = — pAEC. (32) 
Using Eq. ()31|). we get the electric quadrupole moments (J27J) for fractal ellipsoid 

ft 7) = W) + + 7C 2 ) , (33) 



where a = D/3. If D = 3, then we have the well-known relation 

Q(a, p, 7 ) = 9^1 {aA 2 + ^2 + 7C .2 } _ 



VII. CONCLUSION 

In this paper, we have introduced and described the fractional continuous model for the 
fractal distribution of charged particles. Using the fractional calculus and the fractional 
continuous model, we have shown that the fractional integrals can be used for calculation 
of multipole moments of the fractal distribution. The order of fractional integral is equal to 
the fractal dimension of the distribution. 

The fractional continuous models for fractal distribution of particles may have applica- 
tions in plasma physics. This is due in part to the relatively small numbers of parameters that 
define a fractal distribution of great complexity and rich structure. The fractional generaliza- 
tion of integral Maxwell equations may have applications in the analysis of electrodynamical 
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problems involving the fractal stuctures. Therefore, it is interesting to numerically solve the 
fractional equations for charged fractals. The fractional continuous model can be used to 
describe dynamics and kinetics of the fractal distribution in the plasma physics. Extension 
of this model to describe the dynamical properties of fractal distribution by fractional gen- 
eralization of magnetohydrodynamics and Vlasov equations is currently under study by the 
author. 



VIII. APPENDIX 



A. Fractional Gauss's theorem 



In order to realize the representation, we derive the fractional generalization of Gauss's 
theorem 

I (J(r,i),<flS 2 )= I div(J{r,t))dV 3 , 
Jaw Jw 

where the vector J(r, t) = is a field, and div(J) = dJ/dr = dJk/dxk- Here and later 

we mean the sum on the repeated index k from 1 to 3. Using the relation 

2 2 ~ d 

dS d = c 2 (d,r)dS 2 , c 2 (d,r) = T ^ 2 ^ l r l d ~ 2 . 

we get 

/ (J(r,t),dS d )= / C2(d,r)(J(r,*),dS 2 ). 
Jaw Jaw 

Note that we have c 2 (2, r) = 1 for the d = 2. Using the usual Gauss's theorem, we get 

/ c 2 (d,r)(3(r,t),dS 2 ) = [ div(c 2 (d,r)J(r,t))dV 3 . 
Jaw Jw 

The relation 

dV D = c 3 (D, r)dV 3 , c 3 (D, r) = J |r| D ' 3 

in the form dV 3 = 1 (D 1 r)dVD allows us to derive the fractional generalization of Gauss's 
theorem: 

/ (J(r,t),<flS d ) = I c^(D,r)div(c 2 (d,r)J(r,t)) dV D . 
Jaw Jw v ' 

Analogously, we can get the fractional generalization of Stokes's theorem in the form 

/(E,dL r )= / c 2 1 (d,r)(curl(c 1 ( 1 ,r)E),dS d ), 
Jl Js 



where 



2 1 -^T(l/2) , , 
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FIG. 1: Charge of fractal distribution! in the interval [0; 30]. 



B. Cantor set and fractional continuous model 



In the paper, we mention the difference between the real fractal medium structures and 
replacing it by a fractional continuous mathematical model. Some quantitative measure of 
the difference would be helpful. Note that the difference between the real fractal media and 
fractional continuous medium model has an analog of the difference between the real atomic 
structure of the media and the usual continuous medium models of these media. In order to 
have some quantitative measure of the applicability of fractional continuous medium model, 
we can consider the power law for fractal media. The fractal distribution of charged particles 
is characterized by the law Q([0, R}) ~ R D Q- 

The Cantor set is given by taking the interval [0;x], removing the open middle third, 
removing the middle third of each of the two remaining pieces, and continuing this procedure 
ad infinitum. The Cantor set is sometimes also called no middle third set. 

In Fig. 1, we consider the charge y = Q([0,x\) of fractal distribution. The fractal 
distribution is described by the Cantor set with fractal dimension D = ln(2)/ln(3) [22]. 
The continuous model for the charge distribution is described by the continuous line y = x D 
in the interval x G [0, 30]. 

In Fig. 2, we consider the difference between the charge y = Q([0,x]) of fractal distri- 
bution and the charge that is described by fractional continuous distribution in the interval 

x e [0,30]. 
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FIG. 2: Difference between the charge of fractal distribution and fractional continuous model. 
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